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Channel conductance measurements can be used as a tool to study thermally activated electron 
transport in the sub-threshold region of state-of-art FinFETs. Together with theoretical Tight- 
Binding (TB) calculations, this technique can be used to understand the evolution of source-to- 
channel barrier height (E),) and of active channel area (S) with gate bias (V gs ). The quantitative 
difference between experimental and theoretical values that we observe can be attributed to the 
interface traps present in these FinFETs. Therefore, based on the difference between measured 
and calculated values of (i) S and (ii) \dEt/dV gs \ (channel to gate coupling), two new methods 
of interface trap density (Dm) metrology are outlined. These two methods are shown to be very 
consistent and reliable, thereby opening new ways of analyzing in situ state-of-the-art multi-gate 
^ FETs down to the few nm width limit. Furthermore, theoretical investigation of the spatial current 

density reveal volume inversion in thinner FinFETs near the threshold voltage. 
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I. INTRODUCTION 



The non-planar tri-gated FinFET geometry (Fig. [lp,) provides a viable solution to the channel length (L c h) scaling 
due to its better gate to channel electrostatic coupling and reduced Short Channel Effects (SCEs) [JUS]. In a recent 
experimental study of undoped Si n-FinFETs 6J , thermionic emission in the sub-threshold region was used to measure 
(I) the active channel cross-section area (S) (Fig. HI)), which shows the region of channel where the charge prevalently 
flows and (2) the source to channel barrier height (Eb) (Fig.[lp), which reflects on the ease with which electrons travel 
from the source/drain to the channel. Understanding the evolution of the active area 'S' and the barrier height 'E^ 
thus opens up new ways to investigate these FinFETs. 

To shed light into the complicated transport phenomena that can arise in these undoped FinFETs we expand our 
previous work [6j [7] and theoretically investigate the evolution of S and E\>. Since these devices are small and have 
finite number of atoms in the channel, modeling transport requires an atomistic representation of the device. A 20 
band sp 3 d 5 s* atomistic Tight-Binding (TB) model with spin orbit coupling (SO) [BHTU] is well suited for modeling the 
bandstructure of these confined Si channels, since TB can easily take into account the material, geometrical, strain and 
potential fluctuations at the atomic scale [T0|, ITTj . This model also takes into account the coupling of the conduction 
(CB) and the valence bands (VB) which is neglected in simple models like the effective mass approximation (EMA) 
[12] , Thermally activated transport is modeled using a semi-classical 'Top of the barrier' (ToB) model (Fig. [I] c) 
[TUl fT3"] . The simple ToB approach has been shown to model thermionic emission accurately [2] . 

Qualitatively, we found similar theoretical and experimental trends for S versus gate bias (V gs ) and Ef, versus V gs 
[5J. However, the theoretically obtained S and E b values quantitatively over-estimated the experimental values. The 
reduced experimental values can be attributed to the presence of interface traps in these FinFETs [5J [BJ [T51 [To] . 
The effect of interface traps on channel property are even more dominant in the extremely thin FinFETs [7|. This 
difference in S and Eb has been utilized to directly estimate the interface trap density (Du) in these FinFETs, thereby 
eliminating the need to implement special FinFETs geometries to determine D it [T5] . These methods now enable the 
direct implementation of interface trap density metrology in state-of-the-art undoped Si n-FinFETs. 

This paper has been divided into the following sections. Section [TTJ provides the details about the FinFETs for 
which interface trap density metrology has been implemented and the fundamentals of our experimental procedures. 
The details about the self-consistent calculations are provided in Sec. |III A| and the theoretical extraction of Eb and 



S is outlined in Sec. MB Section IV provides the details of the two procedures for obtaining the interface trap 
density. The theoretical and experimental results and the discussion on them are given in Sec. [V] The conclusions 
are summarized in Sec. I VI I 



II. DEVICE AND EXPERIMENTAL DETAILS 



Device details: In this work 7 different FinFETs (labeled A-G) with two different channel orientations of [100] 
((FinFETs A-C and G)) and [110] ((FinFETs D-F)) have been used gj (see Table)!. All the FinFETs have the same 
channel length (L — 40nm). The channel height (H) is either 40nm or 65nm (Table]!]). The channel width (IV) varies 
between 3 to 25nm. All the FinFETs consist of one or more Si nanowire channels etched on a Si intrinsic film with a 
wrap-around gate covering three faces of the channels (Fig. [T]a) [3J. An HfSiO (high-A?) layer isolates a TiN layer from 
the intrinsic Si channel Hj. These FinFETs have either one channel (FinFETs A-C and G) or ten channels (FinFETs 
D-F). These devices have two different surface treatments (with or without H2 annealing) as shown in Table [I] 

Measurement procedure: The experimental value of Eb and S are obtained using a differential conductance (G = 
dlfj/dVds) method. The conductance data are taken at Vd s = V using a lock-in technique. The full experimental 
method and the required ambient conditions have been outlined in detail in Ref. [BJ. 

In the next section details of the theoretical approach to calculate the experimental values of Eb and S in tri-gated 
n-FinFETs are outlined. 



III. MODELING APPROACH 



A. Self-consistent calculation 



The bandstructure for the Si channel is calculated using TB [TUl HU EH- The TB calculation is coupled self- 
consistently to a 2D Poisson solver to obtain the charge and the potential [10J [13] . Once the convergence between the 
charge and the potential is achieved the thermionic current in the FinFETs is obtained using a semi-classical ballistic 
ToB model as shown in Fig. [TJ: [TUJ [131 [T7] . Due to the extensively large cross-section of the devices that combines 
up to 44,192 atoms (for H = 65nm, W — 25nm FET) in the simulation domain, a new NEMO-3D code has been 
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Label 


H 

H 


W 

H 


L 

H 


Channel 
Orientation 


H 2 

anneal 


A 


65 


25 


40 


[100] 


Yes 


B 


65 


25 


40 


[100] 


No 


C 


65 


~5 


40 


[100] 


No 


D 


40 


18 


40 


[110] 


Yes 


E 


40 


18 


40 


[110] 


Yes 


F 


40 


~3-5 


40 


[110] 


Yes 


G 


65 


~7 


40 


[100] 


Yes 



TABLE I. Si n-FinFETs used in this study along with their labels. The surface hydrogen annealing detail is also shown. The 
channel is intrinsic Si, while the source and the drain are n-type doped for all the FinFETs. 



integrated in the top of the barrier analysis [IS]- Since the FinFETs studied here show (i) negligible source-to-drain 
tunneling current and (ii) reduced SCEs [5], the ToB model is applicable to such devices [13]. All the FinFETs are 
n-type doped in the source and drain to a value of 5xl0 19 cm -3 . A 1.5nm SiC>2 cover is assumed. 
Next we outline the procedure to calculate E b and S. 



B. Calculation of Eb and S 



For pure thermionic emission any carrier energetic enough to surmount the barrier from the source (Src) to the 
channel (C) (Fig.JTJ;) will reach the drain (Drn) provided the transport in the channel is close to ballistic |17j . Typically 
Src/Drn in FETs are close to thermal and electrical equilibrium (since heavy scattering in the contacts is assumed 
which leads to instantaneous carrier relaxation) . This allows us to make the assumption that most of the carriers in 
the Src/Drn are thermalized at their respective Fermi-levels (Ef s , Efd in Fig. [l]c). Also the channel potential (U sc f) 
can be determined under the application of V gs using the self-consistent scheme (discussed in Sec III A ) . Hence, for 
the source-to-channel homo-junction inside a FET, the barrier height (Eb) can be determined as a function of V gs , 



E b (V gs ) = U scf (V gs )-E fs . (1) 

This definition of E b implicitly contains the temperature dependence since the simulations are performed at different 
temperatures (T) which enters through the Fermi-Dirac distribution of the Src/Drc. We show in a later section, that 
the temperature dependence of E b in the sub-threshold region is very weak. Therefore, all the theoretical E b results 
shown in this work are at T = 300K. 

The study of thermionic emission model is applicable when the barrier height is much larger than the thermal 
broadening ( E b 3> kgT |19j . where fee is the Boltzmann constant). For this reason, Eq. p\ works only in the 
subthreshold region where E b is well defined [13] • Once the FinFET is above the threshold, E b (< fcgT) is not a well 
defined quantity |13j . Using the E b value, S can be extracted using the conductance (G) in the thermionic emission 
regime for a 3D system [U [19] as, 



G„ = SA-T±ezp(-*jffi) (2) 

where A* is the effective 3D Richardson constant (A* Si3D = 2.1 x 120 Acm~ 2 K~ 2 ) [19], and e is the electronic 
charge. This will hold only when the cross-section size of the FinFET is large enough (i.e.: W, H > 20nm) to be 
considered a 3D bulk system. For a very narrow width FinFET, S cannot be extracted using ^ since the system 
is close to ID. For a ID system the G, under a small drain bias (V^s) at a temperature T, is given by the following 
relation (for a single energy band [2"0]L 



G u 



2e^ 
h 



1 + exp( 



E b (V gs 
knT 



) 



(3) 



where h is the Planck's constant. Since Eq.([3]) lacks any area description, G for ID systems is no more a good 
method to extract S. 
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IV. TRAP EXTRACTION METHODS 



In Ref. [6] it was observed that the active cross-section area (S s im) obtained theoretically was an over-estimation of 
the experimental value (S exp ). In the results section it will be further shown that also the theoretical Ef, value can over 
estimate the experimental Eb value. These mismatches can be attributed to the presence of traps at the oxide-channel 
interface of multi-gate FETs where these traps can enhance the Electrostatic-static screening and suppress the action 
of the gate on the channel [6j [3 [HI [16] . This simple idea is a powerful tool used for the estimation of interface trap 
density (D it ) in these undoped Si n-FinFETs. 



A. Method I: Dm from active area 



Based on the difference between the simulated and the experimental active area (S) values, a method to calculate 
the density of interface trap charges, an, in the FinFETs is outlined. The method is based on the assumption that 
the total charge in the channel at a given V gs must be the same in the experiments and in the simulations. This 
requirement leads to the following, 



S s 



L, 



h ■ Ps 



Sexpt ' L c h ■ p e xpt + e ' &it ' Ech ■ P 



(4) 



where S S i m (S exp t) is the simulated (experimental) active area, P is the perimeter of the channel under the gate (P 
= W + 2H) and p S i m (pexpt) is the simulated (experimental) charge density. By applying Gauss's law at the oxide 
channel interface, p eX pt is obtained from p s %m and an as, 



Pexpt = Psim ~ Pit = Psirn - (e • O lt ■ P) / (W ■ H) 

Using Q and ^ the final expression for an is obtained as, 



(5) 



Vit(V gs ) 



Psim{Vgs)$ sim (Vgs ) 



l - 



e-P 

SexptjVgs) 

(V 38 ) 



n't 



{Vgs) 



WH 



(6) 



[#/cm 2 



This method is useful for wider devices for which Eq.|2]) is valid. For very thin FinFETs (close to a ID system) 
this method cannot be utilized . 

Assumptions in Method I: The extra charge contribution completely stems from the interface trap density (Dit) and 
any contribution from the bulk trap states has been neglected. All the interface traps are assumed to be completely 
filled which implies an = P>n- The interface trap charges are assumed to be situated very close to the oxide-channel 
interface for Eqj5] to be true. Also the interface trap density is assumed to be constant and identical for the top 
and the side walls of the FinFET which is generally not the case [TSJ H] • This method of extraction works best for 
undoped channel since any filling of the impurity/dopant states is neglected in the calculation. Orientation dependent 
Dn for different surfaces could be included as a further refinement. 



B. Method II: D it from barrier control 



The second method does not utilize the Eb value directly but its derivative w.r.t V gs . The term a = \dEb/dV gs \ 
represents the channel to gate coupling [B]. The presence of interface traps weakens this coupling due to the electro- 
static screening. This method of trap extraction is based on the difference in the experimental and the simulated a 
value. The a value can be represented in terms of the channel and the oxide capacitance. The equivalent capacitance 
model for a MOSFET with and without interface traps (Dn) is shown in Fig[2] 

The a value can be associated to the oxide, interface and semiconductor capacitance which is given in Eq.(38) on 
page 383 in Ref. [19]. This leads to the following relation, 

I^I = i-tK W 
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where C+, 



is the total capacitance. For the two cases, as shown in Fig[2] the total capacitance is given by, 

C ox • (Cd + Co) 
Cd + C 
Cd • C ox 



■ ■ 



' Ci, 



'-'tot 



C\i + C 



(8) 
(9) 



where Cu, C ox and Cd are the interface trap capacitance, the oxide capacitance and the semi-conductor capacitance, 
respectively. Eq. ^ represents the capacitance in the experimental device and Eq. ^ represents the capacitance 
in the simulated device under ideal conditions without any interface traps. Combining Eq. Q, Q and ^ and after 
some mathematical manipulations, we obtain, 



Cu 

n ' 

^ or 



Manipulating Eq. (10 1 gives the following relation for Cu, 



Cit — C 



1 



Also Cu can be related to the interface charge density {(Jit) as [T5] . 

d(J it 



Ci, 



dV r 



(10) 



(11) 



(12) 



y* 



In Eq. (Ill all the values are dependent on V gs except C ox - Combining Eq. 11 and 12 and integrating w.r.t V gs 
yields the rmal expression for the integrated interface charge density in these FinFETs as, 



a 



ox 

e 



V2=V T 
VI 

(V gs ) 



( 1 ) 



< exp 



(Vgs) 



- 1 



(Vgs)- 

dV gs [#/cm 2 



where Vt is the threshold voltage of the FinFET and VI is the V gs at which a 
integration range for Eq. ( 13 1 in the sub-threshold region. 



exp/ sim 



(13) 



1. Thus VI and V2 is the 



The second method derived from barrier control has the advantage that it is independent of the dimensionality of 
the FinFET. Hence, Eq. ( 13 1 can be used for wide as well as thin FinFETs. 



Assumptions in Method II: The most important assumption is that the rate of change of the surface potential 
(^(V gs )) is the same as E\, w.r.t V gs . The extra charge contribution completely originates from the density of 
interface trap charges (era) and any contribution from the bulk trap states have been neglected. Also all the interface 
traps are assumed to be completely filled which implies an = Du. This method works best when the change in the 
DC and the AC signal is low enough, such that the interface traps can follow the change in the bias sweep [19] . 



C. Limitations of the methods 



It is important to understand the limitations of the new trap metrology methods to apply them properly. One 
of the main limitation is how closely the simulated FinFET structure resembles the experimental device structure. 
This depends both on the SEM/TEM imaging as well the sophistication of the model. In the present case we 
create the FinFET cross-section structure using the TEM image making the simulated structure as close to the 
experimental device as possible. With the development of better TCAD tools, the proximity of the simulated structure 
to experimental structure has increased. The physical device model needs to comprehend the crystal directions, 
atomistic details, strain and gating realistically to realize the working of the nano-scale FETs. Effective mass models 
fail to properly represent the bandstructure in these type of ultra-scaled nanowires/FinFETs [12j . Our model enables 
the calculation of theoretical conductance value with good confidence to be used in the trap calculation. Furthermore, 
the calculated G is calculated as close to ideal as possible and all the difference between the ideal and experimental 
G is attributed to the traps which may not be true always. An important difference between the two methods is that 
they are calculated over different V gs ranges. This is important since the trap filling and their behavior changes with 
V gs range which should be taken into account accurately. One must also be aware of the embedded assumption of 
complete interface trap filling and the neglect of the bulk traps in the gate dielectric. 
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V. RESULTS AND DISCUSSION 

In this section the theoretical results as well as their comparison with the experimental data are provided. 



A. 3D vs. ID system 

The conduction band structure (E-K) can be utilized to distinguish a ID system from a 3D system. The bulk 
silicon conduction band (CB) has 6 degenerate valleys (Ag) (see inset of Fig. [3] and Fig. |4| which split into 2 sets 
of degenerate valleys called the '4 — 2 configuration' (A 4 — A 2 ) for [100] and [110] ID nanowires (NWs) due to the 
geometrical confinement [TU]. In Si bulk the A 2 valleys are along the [100] direction. For a [100] Si nanowire channel 
the bulk A 2 valleys are projected along the channel axis away from the Y point due to the folding of the Brillouin 
Zone (Fig. [3] a and b) [TU]. The bulk A4 valleys are projected at the Y point [TU]. The bandstructure of conduction 
band for silicon channel with W = 3nm and H = 15nm, and W = H = 15nm is shown in Fig. [3]a and b, respectively. 
For [110] oriented Si channel the valley projection is different compared to the [100] channel. The CB minima is at the 
Off-r position as shown in Fig. [4] This happens due to the different atomic positions and geometrical confinement in 
[100] and [110] channels [TU]. 

The energy separation between the Y and Off-r valleys is given by, 

AE C = E(Y) - E(Off - Y), (14) 

which gives a measure of how close (or far) a ID NW system is from a 3D bulk system. The observation of a large 
AE C value strongly points towards a ID system, whereas a value close to zero points to a bulk system. Tight-binding 
simulations predict a AE C value of around 120 meV for a [100] Si nanowire channel with W = 3nm and H — 15nm 
(Fig{3] a) while this value reduces to 6 meV for a Si nanowires channel with W = H = 15nm (Fig. |3]b). For a [110] 
Si channel the AE C value is around 34 meV for W = 3nm and H — 15nm (Fig|4]a) which reduces to 3.4 meV for W 
= 15nm and H — 15nm (Figfilb). This indicates that larger cross-section silicon channels are closer to the 3D bulk 
system for both [100] and [110] oriented channels. 

The conduction band minimum (CBM) decreases with increasing channel width for both [100] and [110] SiNWs 
(Fig. Fig. [5] a). Also the AE C value decreases with silicon channel width for a fixed height of 15nm ((Fig. [5]b). The 
AE C value is negative for [110] channel since the Off-r valley is lower in energy compared to Y valley. For W > 15nm 
the AE C is less than 5 meV (< k{,T 3 Q x) for both [100] and [110] Si channels. This suggests that silicon channels with 
W > 15nm, and H = 15nm behave electrically close to the bulk Si system at room temperature. 



B. Temperature dependence of Ei 

The source-to-channel barrier height has been assumed to be temperature independent in the sub-threshold region. 
Figure [6] shows the results of a temperature dependent ToB calculations and proves that the barrier height (Eb) is 
only weakly temperature dependent. In the subthreshold region, the Eb value for a device identical to FinFET C, 
is same at four different temperatures ( T=140K, 200K, 240K and 300K). The variation with temperature becomes 
more prominent when the FinFET transitions into the on-state. Thus, Eb has a weak temperature dependence in the 
sub-threshold region allowing us to evaluate Et from the 300K simulations only. 



C. Evolution of E b and S with V ga 

Experimentally, it has been shown that, for undoped silicon n-FinFETs [UJ, Eb reduces as V gs increases. Theoret- 
ically, the Eb value is determined using Eq. |l]) which depends on the self-consistent channel potential (U sc f). As 
the gate bias increases, the channel can support more charge. This is obtained by pushing the channel CB lower in 
energy to be populated more by the source and drain Fermi level [10] . Figure [7] and [8] show the experimental and 
theoretical evolution of Eb in FinFETs G, C and D, E, respectively. Theory provides correct quantitative trend for 
Eb with V gs . Few important observations here are, (i) the theoretical Eb value is always higher than the experimental 
value and (ii) [110] Si devices (D and E) show a larger mismatch to the experimental values. The reason for the first 
point is suggested to be the presence of interface traps in the FinFETs which screen the gate from the channel [6] [7] . 
The second observation can be understood by the fact that [110] channels with (110) sidewalls have higher interface 
trap density due to the higher surface bond density [TU] and poor anisotropic etching of the (110) sidewalls [T51 116] , 
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The active channel area (S) represents the part of the channel where the charge flows [5]. Experimentally S is 
shown to be decreasing with gate bias since the inversion charge moves closer to the interface which electrostatically 
screens the inner part of the channel from the gate [5] . This gives a good indication of how much channel area is used 
for the charge transport. Figure [9] a and b show the experimental evolution of S in FinFET B and E, respectively. 
The theoretical value of S decreases with V gs which is in qualitative agreement to the experimental observation [SJ [7] . 
However, the absolute values do not match. In fact theory over-estimates the experimental S value (Fig. [9]) which is 
attributed to the interface traps. 



D. Trap density evaluation 



In this section we present the results on the interface trap density (Du) in the undoped Si n-FinFETs. 



1. Du using S: Method I 



The calculated D it values for FinFET B and E are 1.02el2cm 2 and 1.81el2cm 2 (Fig 10 a and b, respectively). 
The Du values compare quite well with the experimental Du values presented in Ref.|15j arid also shown in Table 
|llj As expected the D lt value for FinFET E (with [110] channel and (110) sidewalls) is higher than FinFET B ([100] 
channel with (100) sidewalls). This is attributed to the higher D it (~ 2x) on the (110) surfaces (T3JITB]. Our results 
show ~ 1.8x more Du for (110) sidewalls, in close agreement to the experiments. This method allows to obtain the 
Du in the actual FinFETs rather than custom made FETs. 



2. D u using \dE b /dV gs \: Method II 



value needed in this method is taken as ~0.0173 F/m? which is assumed to be the same for all the devices 

The calculated D it values for FinFET C and D are 9.26ellcm -2 

(Fig 



The C ox 

since these FinFETs have similar oxide thickness 

-2 



and 1.563el2cm~ 2 (Figlllla and b, respectively). These calculations also show that the [110] channel device (FinFET 
E) shows a higher Di t compared to the [100] channel device (FinFET C), again consistent to the observations made 
in Ref. [TS] . The advantage of this method is that it can be used to obtain Du in extremely thin FinFETs (close to 



ID system) unlike method I which is applicable only to wider FinFETs (due to the reasons discussed in Sec. MB) 



3. Discussion on the two methods and Du trends 



The D it values for all the FinFETs used in this study are shown in Table [TTJ The important outcomes about the 
two methods arc outlined below: 

• The Du values obtained by the two methods compare very well with the experimental measurement in Ref. |15j 
for similar sized FinFETs (A and B), demonstrating the validity of these new methods. 

• The D it values calculated using method I and II (for B and E) compare very well with each other which shows 
that the two methods are complimentary 7 for large cross-section devices. 

• The Du values calculated for the two similar FinFETs (D and E) compare very well showing the reproducibility 
of the methods [Tj- 

The calculated D it values also reflect some important trends about the FinFET width scaling and surfaces (Table 
|Tlj ). The central points are : 

• Hydrogen passivation considerably reduces Du [5] . This is observed for FinFETs A and B where H2 passivation 
results in - 2x less D lt in FinFET A (7]. 

• Width scaling requires more etching which also increases Du [HI [16] . The same trend is observed in devices A 
to C and D to F (decreasing W). 

• (110) sidewalls show higher Du compared to (100) sidewalls [IS]. The same trend is also observed for FinFETs 
A, B, C, G ((100) sidewall) compared to FinFETs D, E and F ((110) sidewall). 



Device 


Method 


At (10 n cm- 2 ) 


FET type 


Observations 


L=140nm Q2] 


Charge 
Pumping 


1.725 


Special body 
tied FET 




L = 240nm [15] 


2.072 




A 


I 


5.560 


Std. FET 


H 2 
anneal, 
reduces Du 


B 


I 


10.60 


Std. FET 


II 


8.860 


Std. FET 


C 


II 


9.26 


Std. FET 


Thin fin, more Du 


D 


II 


18.31 


Std. FET 


(110) side- wall, 
thin fin, 

more etching, 
more Du 


E 


I 


18.1 


Std. FET 


II 


15.3 


Std. FET 


F 


II 


36.3 


Std. FET 


G 


II 


4.33 


Std. FET 


H2 anneal, less Du 



TABLE II. Values of D it obtained from all the n-FinFETs as well as from Ref.[T5]. 



E. Current distribution 



The spatial current distribution in FinFETs can provide a better insight for optimizing the channel cross-section 
area utilization. Theoretical calculations show that the charge flow in n-FinFETs depends strongly on the geometrical 
confinement. For very small width FinFET (W ~ 5) the entire body gets inverted (Fig[l2]a) and shows little change 
in S with V gs [6]. For wider FinFETs (W ~ 25) the current flow starts from a weak volume inversion and moves 



towards surface inversion as V gs increases (Fig 12 b) [5]. For extremely thin n-FinFETs (W = 5nm, H = 65nm) the 
charge flows through the entire body (volume inversion) compared to the wider n-FinFETs (W — 25nm, H = 65nm) 
where the charge flows at the edges. Thus thin FinFETs can show better channel area utilization for the charge flow 
compared to the wider devices. However thin FinFETs show an increase in Da due to more side-wall etching (Table 
|n| which can severely limit the action of the thin FinFETs. The advancement of fabrication methods and strain 
technology may improve the performance of thin FinFETs as reported by some experimental works [2"Tl - f2"5] . 



VI. CONCLUSIONS 



Two new trap charge density metrology methods in ultra-scaled Si n-FinFETs are presented. The Top-of-thc- 
barrier model, combined with Tight-binding calculations, explains very well the thermally activated sub-threshold 
transport in state-of-the-art Si FinFETs. The qualitative evolution of Ej, and S with V gs are well explained by the 
theory. The systematic mismatch in the experimental and theoretical values of E\, and S led to the development of 
two new interface trap density metrology schemes. The advantage of these schemes is that they do not require any 
special MOSFET structure as needed by the present experimental methods allowing to probe the interface quality of 
the ultimate channel. These methods are shown to provide consistent and reproducible results which compare very 
well with the independent experimental trap measurement results. The calculated trends of interface trap density 
with channel width scaling, channel orientation and hydrogen passivation of the surfaces show a good correlation 
with the experimental observations. Thin width FinFETs could lead to a better channel utilization due to strong 
volume inversion only if surface roughness and the density of interface traps created during the extreme etching of 
these ultra-scaled devices can be reduced. 
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FIG. 1. (a) Scanning-electron-microscope (SEM) image of a Si n-FinFET with [100] channel orientation and single fin. (b) 
The schematic of the cross-sectional cut in the Y-Z plane of a typical tri-gated FinFET. The active cross-section (5*) is in gray, 
H and W are the physical height and width, respectively, (c) Ballistic top of the barrier model employed for calculating the 
thermionic current in the FinFETs. 
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FIG. 2. Equivalent circuits (a) with interface-trap capacitance (Cu) and (b) without interface capacitance. Cd and Cox are the 
depletion and the oxide capacitance, respectively. The idea for this equivalent circuit is obtained from page 381 in Ref. [19] . 
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FIG. 6. Temperature dependence of the simulated barrier height (E b ) in n-FinFET C. At T=300K, V T of the FinFET is 0.62V. 
The overlap of the curves at different temperatures with V gs , below Vr at 300K, shows a weak temperature dependence of Eb 
in the sub-threshold region. The impact of temperature becomes prominent after V gs goes above Vr- 
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